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Abstract. There are few explicit examples in the literature of vector fields exhibiting com- 
plex dynamics that may be proved analytically. We construct explicitly a vector field on the 
three-dimensional sphere S^, whose flow has a spiralling attractor containing the following: 
two hyperbolic equilibria, heteroclinic trajectories connecting them transverselly and a non- 
trivial hyperbolic, invariant and transitive set. The spiralling set unfolds a heteroclinic network 
between two symmetric saddle-foci and contains a sequence of topological horseshoes semicon- 
jugate to full shifts over an alphabet with more and more symbols. The vector field is the 

CO restriction to S"^ of a polynomial vector field in R^. In this article, we also identify global 

^"^ bifurcations that induce chaotic dynamics of different types. 
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^ Any C vector field defined on a compact three-dimensional manifold may be approximated by 

C a system of differential equations whose flow exhibits one of the following phenomena: uniform 

"—I hyperbolicity, a heteroclinic cycle associated to, at least, one equilibrium and/or a homoclinic 
tangency of the invariant manifolds of the periodic solutions — see Arroyo et al [9\. However, 

^ there are no examples of vector fields exhibiting all three features simultaneously 

CO In this article, we construct an explicit example of a C^ vector field on the three dimensional 

^ sphere S^ that is approximated by differential equations exhibiting all the three behaviours. 

ly^ Nearby differential equations may also display heteroclinic tangencies of invariant manifolds of 

^ two equilibria. We describe and characterise some properties of the flow of these differential 

^ equations, whose complex geometry arises from spiralling behaviour induced by the presence of 

Cn saddle-foci. We start with a brief discussion of the literature on examples of this kind. 
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1. Introduction 



h^ 1.1. Lorenz-like attractors. Few explicit examples of vector fields are known, whose flow 

• ^ contains non-hyperbolic invariant sets that are transitive and for which transitivity is robust to 

rS small C^ perturbations. The most famous example is the expanding butterfly proposed by E. 

^ Lorenz (1963) [19j that arises in a flow having an equilibrium at the origin, where the linearisation 

has eigenvalues A^^, Aj, A^ G R\{0} satisfying: 

(1.1) A^ < aJ < < A^ and A^ + AJ > 0. 

In a three dimensional manifold, an equilibrium whose linearisation has real eigenvalues satis- 



fying condition (1.1) is what we call an equilibrium of Lorenz-type. 
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In order to understand the Lorenz differential equations and the phenomenon of robust co- 
existence in the same transitive set of an equihbrium and regular trajectories accumulating on 
it, geometric models have been constructed independently by Afraimovich et al [1] and Guck- 
enheimer and Williams |15j. The construction of these models have been based on properties 
suggested by numerics. 

Morales et al [24j unified the theory of uniformly hyperbolic dynamics and Lorenz-like flows 
stressing that the relevant notion for the general theory of robustly transitive sets is the dom- 
inated splitting. More precisely, they have proved that if A is a robustly transitive attractor 
containing at least one equilibrium of Lorenz-type, then A must be partially hyperbolic with 
volume expanding directions, up to reversion of time. 

In order to discuss this more rigorously, recall that a compact flow-invariant set A is partially 
hyperbolic if there is an invariant splitting TA = E^ ® E^^ for which there are i^, A G R^ such 
that for Vt > 0,Vx G A: 

• \\d,<p{t,x)\Es\\<Ke-'^'; 

• ||a,(/)(t,x)|^.||.||a,(/)(t,x)|^c.^^^J| <Ke-^^ 

where (j){t^p) is the unique solution x{t) of the initial value problem x — /(x), x(0) = p, and 
/ : R^ -^ R^ is a smooth vector field. The direction E^^ of A is volume expanding if 

Mt > 0, Vx G A, det\d:j,(j){t,x)\Egu\ > Ke^\ 

The Lorenz model satisfies these conditions. A good explanation about this subject may be 
found in chapter 3 of Araiijo and Pacifico ^, where it is shown that in a three dimensional 
manifold the only equilibria that exist near a robustly transitive set must be of Lorenz-type ie. 



where condition (1.1) holds - see also Bautista [TT] . 



1.2. Contracting Lorenz Models. In 1981, Arneodo, Coullet and Tresser [7] started the 
study of contracting Lorenz models. The authors considered a variation of the classical Lorenz 



model with respect to the eigenvalues at the origin, in which the condition (1.1) is replaced by: 
(1.2) A^ < AJ < < A^ and A^ + AJ < 0. 

In 1993, Rovella ^29j proved that there exists a contracting Lorenz model x — /(x), x G R^ 
with an attractor A containing an equilibrium so that the following hold: there exists a local 
basin of attraction ;B of A, a neighbourhood U oi f (in the C^- topology) and an open and dense 
subset Uq C U so that for x = g{x) {g G t/o), the maximal invariant set in B consists of the 
equilibria, one or two periodic trajectories, a hyperbolic suspended horseshoe and heteroclinic 
connections. Moreover, Rovella [29] proved that in generic two parameter families x — f(x,/j.), 
with /(^, 0) = /, there is a set of positive measure containing /x = (0, 0) for which an attractor 
in B containing the equilibrium exists. 

The construction of the Rovella attractor [29] is similar to the geometric Lorenz model. Some 
authors constructed contracting Lorenz-like examples through bifurcations from heteroclinic 
cycles and networks: for instance, Afraimovich et al [2] describe a codimension 1 bifurcation 
leading from Morse-Smale flows to Lorenz-like attractors; Morales [23j constructed a singular 
attractor from a hyperbolic flow, through a saddle-node bifurcation. All of these are similar to 



the expanding Lorenz attractor, for which condition (1.1) hold, as opposed to the contracting 
case. 

1.3. Spiralling attractors. Another type of persistent attractor with a much more complex 
geometry occurs near homoclinic and heteroclinic cycles associated to either a saddle-focus or 
a non-trivial periodic solution. This complexity is the reason why the study of this subject was 
left almost untouched for 20 years from L. P. Shilnikov [311 [32] to P. Holmes [16]. Due to the 
existence of complex eigenvalues of the linearisations at the equilibria, the spiral structure of 
the non-wandering set predicted by Arneodo, Coullet and Tresser [8j has been observed in some 
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simulations in the context of the modified Chua's circuit - these attractors are what Aziz-Alaoui 
[lOj cah spiralling attractors. 

Spiralhng attractors were expected for perturbations of homochnic cycles to a saddle-focus 
under a dissipative condition. There is some evidence that in this case spiral attractors exist and 
that they might be persistent in the sense of measure [8J. A sequence of topological horseshoes 
semiconjugate to full shifts on an alphabet with more and more symbols might occur near this 
kind of attractors. 

Nowadays, a particular attention is being given to the study of the dynamics near heteroclinic 
networks with complex behaviour and multispiral attractors. This creates an interest in the 
construction of explicit vector fields whose flows have a specific type of invariant set and for 
which it is possible to give an analytical proof of the properties that guarantee the existence of 
complex behaviour. 

In this article, we describe an explicit example of a polynomial vector field on the three 
dimensional sphere S^ whose flow has a heteroclinic network with two saddle- foci, and a spiralling 
structure containing a hyperbolic non-trivial transitive isolated set. We show that the spiralling 
set is not robustly transitive, but it presents some similarities to those considered by Rovella 

m- 

1.4. Heteroclinic terminology. Throughout this paper, by heteroclinic cycle we mean a set of 
finitely many disjoint hyperbolic equilibria pj (also called nodes)^ j G {1, . . . , A:} and trajectories 
7j, j G {1, . . . , m} such that: 

lim 7-(t) =p^-+i = lim 7j+iW, 

t— )^+oo t— )►— oo 

called heteroclinic trajectories, with the understanding that 7m+i = 7i and Pk+i = Pi- We also 
allow connected n-dimensional manifolds of solutions biasymptotic to the nodes pi and pj, in 
negative and positive time, respectively. In both cases we denote the connection by [pi -^ pj]. 
Throughout this article, the nodes are hyperbolic; the dimension of the local unstable manifold 
of an equilibrium p will be called the Morse index of p. 

A heteroclinic network is a connected set that is the union of heteroclinic cycles, where in 
particular, for any pair of saddles in the network, there is a sequence of heteroclinic connections 
that links them. Heteroclinic cycles arise robustly in differential equations that are equivariant 
under the action of a group of symmetries as a connected component of the group orbit of a 
heteroclinic cycle. We refer the reader to Golubitsky and Stewart [13j for more information 
about heteroclinic cycles and symmetry in differential equations. 

In the next section, we construct an example of a polynomial vector field whose organising 
centre has a heteroclinic network which originates a wide range of phenomena around it. 

2. The vector field and a framew^ork of the article 

Our object of study is the two parameter family of vector fields X(X, Ai, A2) on the unit 
sphere S^ C R"^, defined for X = (a:^i, ^2, X3, X4) E S^ by the differential equation in R^: 

xi = xi(l — r^) — X2 — aiXiX4 + 0^23:^1^4 + A2X3X4 
/p, o\ J ^2 = 3:^2(1 - r'^) + xi - aiX2X4 + a2X2xl 

^ xs = xs{l - r'^) + aiXsX4 + a2Xsxl + X1X1X2X4- X2X1XSX4 

X4 = X4(l — r^) — ai{x'^ — x\ — X2) — Oi2XA{x\ + ^2 + x\) — A1X1X2X3 

where r^ = x^ + X2 + X3 + x\ and a2 < < ai with ai + 0:2 > 0. 



Equation (2.3) with Ai = A2 = was studied by Aguiar et al [3], we describe their results 
briefly. The unit sphere S^ is invariant under the flow of (2.3) and every trajectory with nonzero 



initial condition is asymptotic to it in forward time, so the vector fields X(X, Ai,A2) are well 



SPIRALLING SETS, APRIL 30, 2013 



W^(v) 




X^ = X2 =0 






Figure 1. Dynamics of the organising centre Xo(X) = X(X, 0,0). Left: The one- 
dimensional heterochnic connection from v to w on the invariant circle Fix{SO{2){ji))r\ 
S^. Centre: The two-dimensional heterochnic connection from w to v on the invariant 
two-sphere Fzx(Z2(72)) H S^. Right: Open neighbourhoods of v and w, inside which 
the direction of solutions turning around the connection [v — ^ w] is the same. 



defined. There are two equilibria given by: 

V = (0, 0, 0, +1) and w = (0, 0, 0, -1) 

and the linearisation of X(X, Ai, A2) at (0, 0, 0, e) with 5 = ±1 has eigenvalues 

a2 — SQ^i ^i and a2-\-£Q^i. 

Then, under the conditions above, v and w are hyperbolic saddle-foci, v has one-dimensional 
unstable manifold and two-dimensional stable manifold; w has one-dimensional stable manifold 
and two-dimensional unstable manifold. 

The vector field Xo(X) = X(X, 0, 0) is equivariant under the action of the compact Lie group 
SO(2) e 22(72) where ipo ^ SO(2) acts as: 

'0^(xi, X2,X3,X4) = {Xi COs9 — X2 siuO^Xi Sin0 + X2COS0, X3, X4) 

and 72 G 22(72) acts as: 

72(^1, ^2, Xs, X4) = (Xi, X2, -X3, X4). 

The one-dimensional invariant manifolds of v and w lie inside the invariant circle Fix{SO{2))n 
S^ and the two-dimensional invariant manifolds lie in the invariant two-sphere Fzx(Z2(72)) HS^. 
Thus, symmetry forces the invariant manifolds of v and w to be in a very special position: they 
coincide, see Figure [l] The two saddle- foci, together with their invariant manifolds form a hete- 
rochnic network S that is asymptotically stable by the criteria of Krupa and Melbourne [20|[2T]. 
Indeed since a2 < < ai, it follows that: 



P 



E^ E^ 



a2 — Q^i V 



>1, 



^0^2 + <^i/ 

where Ep and Cp denote the real parts of the expanding and contracting eigenvalues of DXq at 
V and w, respectively. 

The network S can be decomposed into two cycles. Due to the symmetry and to the asym- 
ptotic stability, trajectories whose initial condition starts outside the invariant fixed point sub- 
spaces will approach in positive time one of the cycles. The fixed point hyperplanes prevent 
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5000 



5999 



Figure 2. Left: Projection in the (xa, X4), (xi, X2) and (xi, X4)-planes of the trajectory 
with initial condition (-0.5000,-0.1390,-0.8807,0.3013) for the flow corresponding to 
Xo(X) = X(X, 0,0), with ai = 1, 0^2 = —0.1. Right: Corresponding time series (same 
parameters and same initial condition). The simulations omit the initial transient (the 
independent variable t varies between 5000 and 5999). 



random visits to the two cycles; a trajectory that approaches one of the cycles in S is shown 
in figure [2j The time series of the figure shows the increasing intervals of time spent near the 
equilibria. The sejourn time in any neighbourhood of one of saddle- foci increases geometrically 
with ratio p. 

The symmetries of the organising centre are broken when either Ai or A2 is not zero. In 
sectionjsj after some additional results on Xo(X), we discuss briefly the dynamics of X(X, 0, A2), 
with A2 7^ 0, when the rotation symmetry SO (2) is broken, destroying the network S. Then we 
exploit the more dynamically interesting situation that arises for Ai 7^ and A2 = 0. In this 
case the reflection symmetry Z2 is broken, as well as part of the SO(2)-symmetry. The two- 
dimensional connection breaks into a pair of one-dimensional ones, but due to the remaining 
symmetry, the one-dimensional connection from v to w is preserved. This gives rise to two 
Bykov cycles - heteroclinic cycles with two saddle-foci of different Morse indices, in which the 
one-dimensional invariant manifolds coincide and the two dimensional invariant manifolds have 
an isolated transverse intersection. Analytical proof of transverse intersection of two invariant 
manifolds is usually difficult to obtain but this can be achieved in our example, and is done in 
the Appendix. 

When both Ai and A2 are non zero, all the symmetry is broken, hence the one-dimensional 
connection from v to w disappears. Although the attracting heteroclinic cycles disappears, some 
nearby attracting structures remain, this is discussed in section [4j Our results and conjectures 
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are illustrated by numerical simulations, which have been obtained using the dynamical systems 
packages Dstool and Matlab. 

3. Partial symmetry breaking 

The existence of the attracting heteroclinic network S for Xq requires two separate symme- 
tries to allow structurally stable connections within the two invariant fixed point subspaces. 
Throughout this paper, we control the degree to which the two symmetries are broken by two 
parameters, Ai and A2: Ai controls the magnitude of the SO(2) x Z2(72)-symmetry-breaking 
term and A2 controls the breaking of the pure reflectional symmetry Z2(7i). The parameters Ai 
and A2 play exactly the same role as in [17j. 



Before discussing the effects of breaking the symmetry in (2.3) when either Ai or A2 is non 



zero, we need some additional information on the fully symmetric case. 

3.1. The organising centre. For Xq, there are two possibilities for the geometry of the flow 
around each saddle- focus of the network S, depending on the direction solutions turn around 
[v ^ w], as discussed in Labouriau and Rodrigues [17j . The next proposition shows that each 
solution when close to v turns in the same direction as when close to w. 

Proposition 1. In S^, there are open neighbourhoods V and W of \' and mv, respectively, 
such that, for any trajectory of Xq going from V to W, the direction of its turning around the 
connection [v -^ w] is the same in V and in W. 

Proof. The explicit expression of Xq in spherical coordinates is the special case Ai = of equation 



(A. 6) in the appendix. The equation for the angular coordinate (p in the plane (xi,X2,0,0) is 
(/J = 1. Since this plane is perpendicular to Fzx(SO(2)(7i)), where the connection [v -^ w] lies, 
trajectories must turn around the connection in the same direction. D 

The property in Proposition [l] is persistent under isotopies: if it holds for the organising 
centre Ai = A2 = 0, then it is still valid in smooth one-parameter families containing it, as long 
as there is still a connection. In particular, property (P8) of Labouriau and Rodrigues [17] is 
verified and thus their results may be applied to the present work. 

3.2. Breaking the two-dimensional heteroclinic connection. For Ai 7^ 0, the vector field 
Xi(X) = X(X, Ai,0) is no longer SO(2) © Z2(72) but is still equivariant under the action of 
7i '•— i^TT ^ SO (2). This breaks the two-dimensional connection [w -^ v] into a transverse 
intersection of invariant manifolds. 

Theorem 2. The vector field Xi(X) = X(X, Ai,0) on S^ has symmetry group ^2(71) and for 
small Xi ^ it has a heteroclinic network S'*' involving the two equilibria v and w with the 
following properties: 

(1) there are two one- dimensional heteroclinic connections from vtow inside Fzx(Z2(7i))n 

(2) there are no homoclinic connections to the equilibria; 

(3) the two-dimensional invariant manifolds of v and w intersect transversely along one- 
dimensional connections from w to v. 

Moreover the connections from w to mv persist under small perturbations that preserve the 71- 
symmetry. 

The new network Yj^ is qualitatively different from H. The transversality of VF^(w) and 
VF^(v) and the persistence of the heteroclinic connections from v to w give rise to a Bykov 
cycle. In the non-symmetric context, Bykov cycles arise as bifurcations of codimension 2, the 
points in parameter space where they appear are called T-points in the literature - see figure |3j 
The existence of a Bykov cycle implies the existence of a bigger network: beyond the original 
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Figure 3. Bykov cycle S*: heteroclinic cycle associated to two saddle- foci of different 
Morse indices, in which the one-dimensional invariant manifolds coincide and the two- 
dimensional invariant manifolds intersect transverselly. 

transverse connections, there exists infinitely many subsidiary heteroclinic connections turning 
around the original Bykov cycle. 

Proof. Clearly the perturbation term breaks all the symmetries of Xi(X) except 71. It is also 
immediate that 71 (v) = v, 71 (w) = mv. The circle Fzx(Z2(7i)) H S^ is still flow- invariant and 
contains no other equilibria, thus it still consists of the two equilibria v and mv and the two 
connections from v to mv. The plane Fzx(Z2(7i)) will remain invariant under any perturbation 
that preserves the symmetry. Since on this plane v is a saddle and w^ is a sink, the connection 
persists under small perturbations. Since v and w are hyperbolic, Fzx(Z2(7i)) D [v ^ w], and 
dimVF^(v) = dimVF^(w) = 1, item S follows. Breaking the Z2(72)-equivariance is necessary 
for the existence of transverse intersection of the manifolds VF^(w) and VF^(v). Nevertheless, 
the manifolds could intersect non transversely. The proof of transversality of the intersection of 
the two-dimensional manifolds, using the Melnikov method is similar to that given in Aguiar et 



al [3j. For completeness, we give the proof in Appendix (A.2). D 



When Al ^ 0, we observe an explosion of suspended horseshoes in the neighbourhood of U'*', 
a phenomenon often called instant chaos. 

Theorem 3. For generic ai ^ 0:2 and for any small Ai 7^ 0; the vector field Xi = X(X, Ai, 0) 
is 7i2{")i)-equivariant and has a compact spiralling set A C S^ containing the following: 

(1) a heteroclinic network Y,'^ involving the two saddle- foci v and w with a transverse inter- 
section of the two dimensional invariant manifolds; 

(2) a uniformly hyperbolic compact set % — [j-^^^Hi, where {Hi}i^z is an increasing se- 
quence of invariant sets, accumulating on Yj^ . Each Hi is topologically conjugate to the 
suspension of a full shift over a finite set of symbols. 

Proof. The proof of ([l]) follows from Theorem [2J Note that transversality implies the existence 
of a larger network with infinitely many heteroclinic connections. 

For Al = 0, the heteroclinic network S for Xq is asymptotically stable. In particular, there 
are arbitrarily small compact neighbourhoods of S such that the vector field is transverse to 
their boundaries where it points inwards. Let M be one of these neighbourhoods. There exists 
A^ > such that for Ai < A^, the vector field Xi is still transverse to the boundary of AA and 
S"^ C N . Thus AA is a compact set that is positively invariant under the flow of Xi. Hence A/" 
contains an attractor. 

If Q^i ^ 0:2, using Samovol [30^ the flow of Xi is C^-linearisable around each saddle- focus. 
Linearisation may fail under resonance conditions that correspond to curves in the («i,«2)~ 
plane. This restriction has zero Lebesgue measure, and thus it does not place serious constraint 
on our analysis. We use the basic ideas of Aguiar et al [4J ^iid Aguiar et al [Sj adapted to our 



SPIRALLING SETS, APRIL 30, 2013 




Figure 4. Construction of the suspended horseshoe Hi for Xi near the cycle, giving 
rise to the Cantor set on the wah of the cyhnder, generated by an arbitrarily large number 
of strips. 



purposes to show the existence of the suspended horseshoe T-L accumulating on the heteroclinic 
network (see figure B. Inside each open set where the linearisation has been done, we take a 
cylindrical neighbourhood of each equilibrium, V and W. 

The boundary of each cylinder forms an isolating block: the flow is transverse to the cylinder 
walls, top and bottom. Near the saddle-focus v, the flow goes in at the cylinder walls, and it 
goes out at the top and bottom. Near w, the flow goes in at the cylinder top and bottom and it 
goes out at the wall. Inside the cylinder the vector field is linear, so the transition from the wall 
to top/bottom and from top/bottom to the wall is well understood. The transition between 
cylinders takes place in a flow-box fashion around each connection. 

Any segment of initial conditions lying across the stable manifold of v is wrapped by the flow 
around the isolating block. By the A-lemma it accumulates as a spiral on the unstable manifold 
of V and then on the stable manifold of the next saddle, w - see figure |4j This spiral of initial 
conditions on the top/bottom of the neighbourhood around w is mapped, by the local map near 
w, into points lying on a helix accumulating on its unstable manifold, which crosses transversely 
the local stable manifold of v infinitely many times. 
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The composition of consecutive local maps and transition functions, transforms the original 
segment of initial conditions lying across VF^(v) into infinitely many segments with the same 
property. Thus, it is possible to construct recursively a nested sequence of strips accumulating 
on the stable manifold of v that follow prescribed sequences of heteroclinic connections and that 
return to the cylindrical neighbourhood of v. Using now Proposition 3 and Corollary 5 of Aguiar 
et al [5 adapted to our purposes, instead of a nested sequence of intervals, we may construct, 
as in figure [4| a sequence of rectangles 7^^, accumulating on S'*', such that the forward iterate of 
initial conditions in TZi lie in strips Si that cross transversely all the TZi, giving rise to a chain 
of suspended horseshoes Hi. D 

The suspended horseshoe T-L containing the heteroclinic network on its closure is semi- conjugate 
to a full shift over an infinite alphabet with countably many symbols. The overall dynamics is 
organised around the suspended horseshoes. In [3], it is shown that that if 11 is a transverse 
section to the suspended horseshoe, the first return map of 7{ to 11 is hyperbolic at all points 
where the horseshoe is well defined. Points lying on the invariant manifolds of v and w are 
dense in T-L. This chain of horseshoes might have positive Lebesgue measure like the "fat Bowen 
horseshoe". Rodrigues [27j proved that the shift dynamics does not trap most trajectories in 
any small neighbourhood of S'*': 

Corollary 4. Let N^ be a tubular neighbourhood of one of the Bykov cycles S'*'. Then in any 
cross-section Sq C N^ at a point q in\w ^ v] the set of initial conditions in Sq H A^^ that do 
not leave N^ for all time has zero Lebesgue measure. 

3.2.1. Heteroclinic Switching and Subsidiary Dynamics. The overall dynamics is organised around 
the suspended horseshoes. Given a cross section E to S'*', in the restriction to an uniformly hy- 
perbolic invariant compact set, the dynamics can be conjugated to a full shift over a finite 
alphabet. In particular, the topological entropy of the corresponding flow is positive. One aston- 
ishing property is the possibility of shadowing the heteroclinic network S'*' by the property called 
heteroclinic switching: any infinite sequence of pseudo-orbits defined by admissible heteroclinic 
connections can be shadowed, as we proceed to define. 

For the heteroclinic network Ti'^ with node set {v,vv^}, a path of order k on S'*' is a finite 
sequence s^ = {(^j)je{i,...,k} of heteroclinic connections Cj = [Aj -^ Bj] in T^~^ such that Aj, Bj E 
{v,vv^} and Bj = ^j+i- An infinite path corresponds to an infinite sequence of connections in 
S'*'. Let A^j]* be a neighbourhood of the network S'*' and let Ua be a neighbourhood of each node 
A. For each heteroclinic connection in S'*', consider a point p on it and a small neighbourhood 
V oi p. We assume that the neighbourhoods of the nodes are pairwise disjoint, as well for those 
of points in connections. 

Given neighbourhoods as before, the trajectory (pf{t, g), follows the finite path s^ = {cj)je{i,...,k} 
of order A:, if there exist two monotonically increasing sequences of times (t^)^^{i /^.^i} and 
{^i)ie{i,...,k} such that for alH G {1, . . . , k}, we have ti < Zi < ti^i and: 

(1) (ff{t,q) C A^E- for ah t E (ti,t/e+i); 

(2) ipf{ti,q) G Ua^ and (p{zi,q) G Vi and 

(3) for all t ^ {zi^ ^i+i), ^/(^5 q) does not visit the neighbourhood of any other node except 
that of Ai^i. 

There is finite switching near S"*" if for each finite path there is a trajectory that follows it. 
Analogously, we define infinite switching near S'*' if every forward infinite sequence of connections 
in the network is shadowed by nearby trajectories. Infinite switching near S'*' follows from the 
proof of Theorem Is] and the results of Aguiar et al [5j . The solutions that realise switching lie 
in a tubular neighbourhood A^^ of S'*', hence from the results of Rodrigues |2Z|, it follows that: 
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Figure 5. The spiralling set A. Top: Projection in the (x3,X4), (xi,X3) and (xi^x^)- 
planes of the trajectory with initial condition (-0.5000, -0.1390, -0.8807, 0.3013) for the 
flow of Xi, with ai = 1, Qf2 = —0.1 and Ai = 0.05. Centre and Bottom: Projection 
in the (x2,X3,X4), (xi,X2,X3), (xi,X2,X4) and (xi, X3, X4)-hyperplanes of the trajectory 
with the same initial condition and parameters. 



Corollary 5. There is infinite switching near the network T^'^ . Infinite switching is realized by 
a set of initial conditions with zero Lebesgue measure, whereas finite switching occurs for a set 
of infinite conditions with positive Lebesgue measure. 



The complex eigenvalues force the spreading of solutions around the unstable manifold of w, 
allowing visits to all possible connections starting at w. The transversality enables the existence 
of solutions that follow heteroclinic connections on the two different connected components of 
S^\VI/^^^(v), the upper and the lower part on the wall of the cylinder - see figure ffl 

From the results of [TT], it also follows that near S'*' the only heteroclinic connections from v 
to w are the original ones, that there are no homoclinic connections associated to the equilibria 
and that the finite switching near S'*' may be realised by an n-pulse heteroclinic connection 
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Figure 6. Time series for the trajectory with 
(-0.5000,-0.1390,-0.8807,0.3013) for the flow of Xi, with ai 
and Ai = 0.05. 



initial condition 
= 1, a2 = -0.1 



[w ^ v]. Moreover, among all the solutions which appear in T-L^ there are infinitely many knot 
types, inducing all link types. 

We illustrate the chaotic behaviour in the projeted phase portraits of figure [5] and in the 
time series of figure [6] corresponding to a trajectory that stays near the heteroclinic network. 
Observing figure [5J it is clear why we call the nonwandering set A of Xi a spiralling set. The 
time series of figure [6] also suggest heteroclinic switching: the trajectory follows a sequence of 
heteroclinic connections in a random order. We stress that this can occur because the symmetries 
have been broken. 

Under generic perturbations (not necessarily equivariant), there is still an invariant topological 
sphere, since S^ is normally hyperbolic, and the transverse connections are preserved. The 
spiralling set presented in Theorem [3] may lose branches, changing its nature. Any compact 
flow-invariant set near the ghost of the heteroclinic cycle, disjoint from it, contains an uniformly 
hyperbolic basic set displaying structurally stable homoclinic classes in its unfolding, where the 
results of Rodrigues et al [28j may be applied. We will return to this issue in section [4] below. 

3.3. Nonhyperbolic behaviour. Based in [iHj, applying the construction of the proof of The- 
orem [3] to the unstable manifold of w, we obtain the following result: 

Theorem 6. There are values X\ arbitrarily close to 0, for which the flow of^(X, X\, 0) has a 
heteroclinic tangency between VF^(w) and W^{'v), coexisting with the transverse connections in 

The tangencies of Theorem |6] take place far from the invariant set T-L of Theorem [3| Newhouse 
results [25j 26j ensure the existence of infinitely many sinks nearby. 



Proof. First of all, we define some terminology. As in the proof of Theorem pi (see figure 10) 



let V and W be cylindrical neighbourhoods of each equilibrium point after linearising X around 
them. Let 

• O^^^ be the top/bottom of the cylinder V, where the flow goes out of V; 

• I^ be the wall of the cylinder V^ where the flow goes in V., 
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Figure 7. Left: Projection in the (xi, X2), (^2, xs) and (^3, X4)-planes of the trajectory 
with initial condition (-0.5000,-0.1390,-0.8807,0.3013) for the flow corresponding to 
the equation (2.3), with Ai = 0, A2 = 0.05 for ai = 1 and 0^2 = —0.1. The dotted line 
on the (x3,X4)-plane indicates the position of the original cycle. Right: Time series for 
the corresponding trajectory (same initial condition and same parameters). 



r\out 
Tin 



be the wall of the cylinder W ^ where the flow goes out of VK; 
be the top/bottom of the cylinder VF, where the flow goes in W . 



Consider the closed curve in where WJ^^^iyv) intersects In('v), as shown in figure 10, For small 
Ai 7^ 0, this curve has a point of maximum height that divides it in two components, each one 
of which is mapped into a helix around Out{'w). The two helices taken together form a curve 
with a fold point (see Figure [s]). Varying Ai moves the fold point around Oixt(w), exponentially 
fast in Ai. As Ai tends to zero, the fold point will cross Wf^^{'v) nOut(-w) infinitely many times, 
creating heteroclinic tangencies. 

D 



3.4. Breaking the one-dimensional heteroclinic connection. For A2 7^ 0, the vector field 
X2 = X(X, 0, A2) is no longer SO(2)-equi variant but is still equivariant under the action of 
72- Since the heteroclinic connections [v -^ w] lie on Fzx(Z2(7i)) = {(0, 0, X3, X4), X3, X4 G R}, 
these connections disappear, destroying the heteroclinic network S, but for small values of A2 
there will still be an attractor lying close to the original cycle. Theorem 6 of [17j shows that 
for sufficienly small A2 7^ 0, each heteroclinic cycle that occurred in the fully symmetric case 
is replaced by a stable and hyperbolic periodic solution. Using the reflection symmetry 72, 
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Out(w) 




Figure 8. Heteroclinic tangency between W^(w) and W*(v), coexisting with the 
transverse connections in S"^. 



w-homoclinics 

/'I 





Figure 9. Left: Quahtative bifurcation diagram for the family X(X, Ai,A2). At the 
sohd curves tangent to the Ai axis there is an expanding Shilnikov homochnic at w. 
Consecutive pairs of these curves bound a tongue (shaded) where X has infinitely many 
unstable periodic solutions associated to a suspended horseshoe. Contracting Shilnikov 
homoclinics at v appear at the dotted curves tangent to the A2 axis, that limit tongues 
where X has an attracting period orbit. Right: Two linked Shilnikov homoclinic loops 
appear when the two types of curve meet. 

two stable periodic solutions co-exist, one in each connected component of S^\Fzx(Z2(72)). 
Their period tends to +oc when A2 ^ and their basin of attraction must contain the basin of 
attraction of S. For sufficiently large A2, saddle- node bifurcations of these two periodic solutions 
occur. Figure [7| illustrates the existence of a single attracting periodic solution in each connected 
component S^\Fzx(Z2(72)) of the phase space. 



4. Breaking all the symmetry 

4.1. Linked homoclinic cycles. In this section, we describe the global bifurcations that are 
most important for our analysis, which appear when both Ai and A2 are non-zero and X has 
no symmetry. Different chaotic dynamics in (2.3) organise a complex network of bifurcations 
involving rotating nodes, that had not yet been considered. The next theorem describes the 
bifurcation diagram for X, shown in Figure [9j 

Theorem 7. The germ at the origin of the bifurcation diagram on the (Xi^ X2) -plane for 
X(X, A1A2) when A1A2 7^ satisfies: 
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(1) For each (Ai, A2); there is a suspended uniformly hyperbolic horseshoe, topologically con- 
jugate to a full shift over a finite number of symbols. As (Ai, A2) tends to (Ai,0) with 
Ai 7^ 0^ the horseshoe accumulates on the set % of Theorem\^ 

(2) There exists a countable family of open tongues, tangent at (0, 0) to the \2-axis, where X 
has an attracting periodic orbit. The closures of these tongues are pairwise disjoint. At 
the two curves comprising the boundary of each tongue there is an attracting Shilnikov 
homoclinic connection at v. 

(3) There exists a countable family of open tongues, tangent at (0, 0) to the Xi-axis, where X 
exhibits a suspended nonhyperbolic horseshoe, topologically conjugate to a full shift over 
an finite set of symbols. The closures of these tongues are pairwise disjoint. At the two 
curves comprising the boundary of each tongue there is an unstable Shilnikov homoclinic 
connection at w. 

(4) The two sets of curves in the parameter plane described in ^ and ^ meet at isolated 
points where the two types of homoclinic connections coexist. These codimension two 
points accumulate at the origin. The closure of the pair of homoclinic orbits forms a link 
whose linking number tends to infinity as the points approach the origin of the (Ai, A2)- 
plane. 

Theorem [t] follows from Theorems 7 and 8 in [17} . We sketch here an alternative proof, 
emphasising the transition between the two different types of chaos and their geometry. 

Proof We use the sets 4^, 0^^^ /^ and 0^\ defined in the proof of Theorem g In Theorem g 
for A2 = 0, we have found a nested chain of suspended horseshoes H contained in the spiralling 
set A. When generic perturbation terms are added, the hyperbolic horseshoes T-L lose infinitely 
many branches, only a finite number of strips survive, as stated in assertion ([l]). When Ai 7^ 
and A2 = 0, a segment s of initial conditions in /^^, with s transverse to VF^(v), is mapped into a 
spiral in I^ and then into a helix in O^ ^ accumulating on VF^(w). When A2 7^ 0, the spirals in 
I^ are off-centred (case (a) of figure [ll|) and will turn a finite number of times around W^{w). 
This is why the suspended horseshoes in assertion ^ have a finite number of strips. Suspended 
horseshoes arise when VF^(v) is close to VF^(w) H /^, the number of branches increases when 
VF^(v) approaches VF^(w) H/^. 

From the coincidence of the invariant two dimensional manifolds near the equilibria in the 
fully symmetric case, we expect that when Ai is close to zero, VF^(v) intersects the wall O^^ of 



the cylinder VF in a closed curve as in figure 10 When both symmetries are broken, the Bykov 



cycle Ti^ is destroyed, giving rise to Shilnikov homoclinic cycles involving the saddle-foci v and 
w. The equilibria have different Morse indices, thus the dynamics near each homoclinic cycle 
is qualitatively different, see Shilnikov [31l [32]: all the homoclinic orbits associated to v are 
attracting since the contracting eigenvalue is larger than the expanding one; Each homoclinic 
cycle associated to w has a suspended horseshoe near it and thus infinitely many periodic 
solutions of saddle type occur in every neighbourhood of the homoclinicity of w. 

The existence of homoclinic orbits is not a robust property, they occur along the curves in 
the (Ai, A2) plane described in assertions ^ and ^ - details on these curves are given in [T7] . 
At values of (Ai, A2) where the two types of curves cross, as in assertion Q, the homoclinics at 
V and w occur at different points in phase space. 



Case (b) of figure 11 describes the homoclinic cycle at w: initially VF^(w) meets I^ at a 
closed curve, that is then mapped into a double spiral in I^ centred at VF^(v) H/^. When this 
spiral meets VF^(w) H /^, an unstable homoclinic cycle associated to w is created. Similarly, 
the backwards image of VK^(v) in O^^ is a closed curve, that iterates back into a double spiral 
in I^ centred at VF^(w) H/^, creating an asymptotically stable homoclinic orbit of v when this 
spiral meets VF^(v) H I^. The values of (Ai, A2) where these intersections happen correspond 
to the two types of curves in figure^ The double spiral VF^(w) H I^ is mapped into a spiral in 
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V 



w 



Figure 10. Conventions in the cylindrical neighbourhoods of v and w used in Theo- 
rem m Both Wf^^{v) n O^^ and Wl^^^{-w) D OIP are closed curves in the boundaries of 
the cylinders for small values of Ai ^ 0. 



O^^ as in figure 11 (c). As one of the spiral arms gets closer to W^{w) in /^, its image winds 
around 0^^ 

Along the curve where one of the arms of the double spiral in VF^(w) H I^ touches VF^(w) 
in /^, there are some parameter values where the centre of the spiral in O^^ meets the closed 
curve VF^(v) H O^^ (figure 
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(c)). At these points both homoclinic cycles coexist in different 
regions of the phase space and the basin of attraction of the homoclinic at v lies inside the gaps 
of the suspended horseshoes that appear and disappear from the homoclinicity at w. D 

We observe a mixing of regular and chaotic dynamics, in regions of parameter space (Ai, A2) 
where the two types of tongues overlap. Examples of chaotic behaviour of this type should be 
regarded as quasiattractors due to the presence of stable periodic orbits within them. A lot 
more needs to be done before these transitions are well understood. 

4.2. Finite heteroclinic switching. Generic breaking of the 71-symmetry destroys the het- 
eroclinic network, because the two connections [v -^ w] break. Nevetheless, finite switching 
might be observed: for small values of A2 there are trajectories that visit neighbourhoods of 
finite sequences of nodes. This is because the spirals on top of the cylinder around w are off- 
centred and will turn a finite number of times around VF^(w). In this way we may obtain points 
whose trajectories follow short finite paths on the network. As VF^(v) gets closer to VF^(w) (as 
the system moves closer to symmetry) the paths that can be shadowed get longer. 

5. Discussion 

The goal of this paper is to construct explicitly a two-parameter polynomial differential equa- 
tion, in which each parameter controls a type of symmetry breaking. We also prove analytically 
the transverse intersection of two invariant manifolds using the adapted Melnikov method which 
is very difficult in general. Some dynamical properties follow applying results in the literature. 
We also describe the transition between the different types of dynamics. As far as we know, our 
construction is new. 

For three-dimensional flows, one homoclinic cycle of Shilnikov type is enough to predict the 
existence of a countable infinity of suspended spiralling horseshoes. Their existence does not 
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Case (a) 



Case (b) 




W (w) 



W"(v). 




W (w) 





Case (c) 




upper part of Q^^ 




-^^^ 


w» 




^^^ 


< W^(w) 


W'(v) 




V] \ 





W'(w) 



[v ►w] 





Wlv) 



Figure 11. There is a sequence of values of A2 for which horseshoes with infinitely 
many legs appear and disappear as in a dance, being part of a very complex spiralling 
network of sets near the ghost of E*. (a) A segment in I^ transverse to W^{v) is 
mapped into a spiral in I^. When A2 varies, branches of this curve get closer to W^{w) 
and its image makes more turns around O^^, creating horseshoes with more strips, (b) 
When one of the arms of the double spiral W'^{'w) H I^ touches VK^(w) a homoclinic 
loop at w is created, (c) The double spiral VK^(w) n I^ is mapped into a curve in O^^ 
that winds a few times around the cylinder and then spirals into a point in VK^(v). A 
stable homoclinic loop at v is created when the spiral centre touches the closed curve 
W^{v)r]0^*. This may happen for the same parameter values that yield the homoclinic 
at w. 



require breaking the homoclinic connection as in the case of saddles whose linearisation has only 
real eigenvalues. The study of spiralling sets has not attracted as much attention as the Lorenz 
attractor because, in general, it is difficult to understand the topology associated to non-real 
eigenvalues. 

The spiralling set A of Theorem [3] is remarkably different from the robustly transitive sets 
described by Araiijo and Pacifico [6J and Morales el al [24j. This is easily seen using assertion ([s]) 
of Theorem [7| to obtain, for arbitrarily small A2, a homoclinic cycle to w of Shilnikov type 
with the expanding eigenvalue greater than the contracting one. This induces the appearance 
of an arbitrarily large number of attracting or repelling periodic solutions. If A were robustly 
transitive, this would imply that Xi could be C^-approximated by vector fields exhibiting either 
attracting or repelling sets, which would be a contradiction. 

Although the spiralling set is not robustly transitive, by Theorem [2] it is persistent under 
Z2-symmetric perturbations. The fact that these wild spiral sets contain equilibria makes them 
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similar to Lorenz-like attractors described by Rovella [29j. Both examples show suspended 
horseshoes, coexisting with equilibria and periodic solutions. In our example, suspended horse- 
shoes are present for an open set of parameters whereas in Rovella's example they appear for 
a set of parameters with positive Lebesgue measure. Moreover, in our example, the suspended 
horseshoes coexist with heteroclinic tangencies and these in turn give rise to attracting periodic 
orbits that may be situated in the gaps of the horseshoes. The attractor splits into infinitely 
many components, whereas Rovella's example satisfies Axiom A for a set of parameters with 
positive Lebesgue measure. 

A generic perturbation of the partially symmetric vector field Xi still has some spiralling 
structure. It contains a finite sequence of topological horseshoes semiconjugate to full shifts 
over an alphabet with a finite number of symbols, instead of an infinite sequence. The dynamics 
near these specific networks can be lifted to larger networks. 

The properties of the family X that we have described depend strongly on the orientation 
around the heteroclinic connection [v ^ vvr]. If the rotations inside V and inside W had opposite 
orientations, for some trajectories these two rotations in V and W would cancel out. Newhouse 
phenomena would occur near a Bykov cycle of this type, and hence infinitely many attractors. 
A systematic study of this case is in preparation. 

A lot more needs to be done before we understand well the dynamics of systems close to 
symmetry involving either saddle-foci or rotating nodes. Besides the interest of the study of 
the dynamics arising in generic unfoldings of an attracting heteroclinic network, its analysis 
is important because in the fully non-equivariant case the explicit analysis of the first return 
map seems intractable. In this article we were able to predict qualitative features because the 
non-symmetric dynamics is close to symmetry. 

Acknowledgements: The authors would like to express their gratitude to Manuela Aguiar 
(University of Porto) for helpful discussions at the beginning of this work. Also special thanks 
to Maria Luisa Castro (University of Porto) for the numerical simulations in Matlab shown in 
figure |5) 

Appendix A. Transversality of invariant manifolds 

A.l. Melnikov method revisited. Melnikov [22^ studied a method to find the transverse 
intersection of the the invariant manifolds for a time periodic perturbation of a homoclinic 
cycle. The pioneer idea of Melnikov is to make use of the globally computable solutions of 
the unperturbed system the computation of perturbed solutions. We start this appendix with 
a short description of this theory applied to saddle-connections. For a detailed proof for the 
homoclinic case, see Guckenheimer and Holmes |14^ section 4.5]. 
If X G R^, t G R and < 6 « 1, consider the planar system: 

(A.4) X = fiX) + sg{X,t) 

such that: 

• there exists T > such that g{X^ t) = g{X, t + T) for all t ie g is T-periodic; 

• for 6 = 0, the flow has a heteroclinic connection Fq associated to two hyperbolic points 
po andpi; 

• the unstable manifold of po coincides with the stable manifold oi pi. 
Associated to the system (A.4), taking S^ = R/^, we define the suspended system: 

X = f{X) + eg{X, 9), ^ = 1, (X, 0) G R^ x S^. 

With the above assumptions, {po} x S^ and {pi} x S^ are hyperbolic periodic solutions for 
the suspended flow, whose invariant manifolds Fq x S^ coincide. Since the limit cycles are 
hyperbolic, their hyperbolic continuation is well defined for 6: ^ 0; hereafter we denote them 
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by {pq} X S^ and {pf} x S^. Under general conditions (without symmetry for example), the 
heteroclinic connection Fq x S^ is not preserved and for a non-empty open set in the parameter 
space, the invariant manifolds meet transversely. Their intersection consists of a finite number 
of trajectories. 

The splitting of the stable and unstable manifolds on a transverse cross section S for the 
perturbed system is measured by the Melnikov function: 

(A.5) M(to) = f f{qo{t)) A g{qo{t),t + to), exp (- f trDf{qo{s))dsj dt, 

where qo{t) is the parametrisation of the solution of the unperturbed system ( A.4| ) starting at 
to = on To . Recall that the wedge product in R^ of two vectors {ui, U2) and (i^i, i;2) is simply 
given by U1V2 — U2V1. When / is hamiltonian, then trDf(qo(t)) = and the Melnikov function 
is simpler. The main result we will use is the following: 

Theorem 8 (Bertozzi ^2j, Melnikov ^, adapted). Under the above conditions, and for e > 
sufficiently small, if M{to) has simple zeros, then W^(pq) and W^{pi) intersect transversely. 

This result is important because it allows to prove the existence of transverse (and hence 
robust) homo and heteroclinic connections as we proceed to do. 

A. 2. Proof of transversality. 

Theorem 9. //Ai 7^ and A2 = 0^ the two dimensional invariant manifolds of the equilibria w 
and V o/X(X, A1A2) intersect transversely in S^ along one- dimensional orbits. 

Proof. In spherical coordinates 

xi = r sin (j) sin 9 cos (f X2 — r sin (j) sin 6 sin if x-^ — r cos (j) sin 6 X4^ — r cos 6 



equations (2.3) restricted to S*^ can be written as: 

9^ai sin6>cos(2(/)) + ^ sin(26>) + ^ s\n^{(j)) s\n^{9) cos((/)) sin(2(/^) 
(A.6) = -ai cos(6>) sin(20) - ^ sin3(0) sin(2(9) sin(2(/^) 



From the equation (/3 = 1, we get Lp{t) —t^t^Yi and (A.6) is reduced to a three-dimensional 
non- autonomous differential equation of the form: 

<i>^f2{0,(t^) + Xi92{eA,t) 

where 

/i(6>,(/)) = aisin6>cos(2(/)) + ^ sin(26>) /2(6>,0) = -ai cos(6>) sin(2(/)) 

gi{e,(t),t) = ^sin2((/))sin2(6>)cos((/))sin(2t) g2{0,(t),t) = -| sin^((/)) sin(26>) sin(2t) . 

The maps gi and ^2 are periodic in t of period tt. For the Melnikov function M{to) defined 
in (A.5) we write / = (/i,/2) and g — (^1,^2)- The parametrisation qo{t) of the connection 
[w -^ v] in the unperturbed system, Ai = 0, is defined by = | + /ctt, k G {0, 1}. Thus, in the 
unperturbed system, the connections [w -^ v] are parametrised by: 

Therefore, for A; G {0, 1}, we have: 

fMit)) = (-l)'=aisin(^(t)) + f sm(2^(t)) h{qi{t)) = 

gi{ql{t),t + to) = g2{ql{t),t + to) = (-l)'=+4 sm(2^(t)) sm(2(t + to)) • 
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To see that the integral M{to) converges, note that the exterior product in the definition of 
the Melnikov function M{to) is bounded since it is given by: 

= [-ai sm{9{t)) + (-1)^+1 f sin(20(t))] | sm{29{t)) sin(2(t + to)) 

and 

trDf{qo{s)) = ai cos{e{s)) + 0^2 cos(26>(5)) 
hence the Melnikov integral M{to) does not depend on Ai. Since it has been shown in Aguiar 
et al ^ lemma 16] that for any r > 0, the integral 

/ exp f — / a2r cos{9{s)) + a'^r'^ cos{29{s))ds \dt 

converges, the convergence of M(to) follows. Finally from the proof of Proposition 12 in Aguiar 
et al [3] it follows that the roots of M(to) are simple, completing the proof of theorem [9J D 
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